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Abstract 

Water wave propagation can be attenuated by various physical mechanisms. One 
of the main sources of wave energy dissipation hes in boundary layers. The present 
work is entirely devoted to thorough analysis of the dispersion relation of the novel 
visco-potential formulation. Namely, in this study we relax all assumptions of the 
weak dependence of the wave frequency on time. As a result, we have to deal with 
complex integro-differential equations that describe transient behaviour of the phase 
and group velocities. Using numerical computations, we show several snapshots of 
these important quantities at different times as functions of the wave number. Good 
qualitative agreement with previous study [Dut09] is obtained. Thus, we validate 
in some sense approximations made anteriorly. There is an unexpected conclusion 
of this study. According to our computations, the bottom boundary layer creates 
disintegrating modes in the group velocity. In the same time, the imaginary part of 
the phase velocity remains negative for all times. This result can be interpreted as 
a new kind of instability which is induced by the bottom boundary layer effect. 

Key words: free-surface flows, visosity, dissipation, visco-potential flow, bottom 
boundary layer 



1 Introduction 



The classical potential free-surface flow theory is known to be a good and rel- 
atively inexpensive model of water waves (especially in comparison with free- 
surface Navier-Stokes equations formulation [HW65, SZ99, WY07]). However, 
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there are some physical situations where viscous effects cannot be neglected. 
The necessity of including some dissipation into various water waves models 
was pointed out explicitly in a number of experimental studies [ZG71, Wu81, 
BPS81]. For example, in the "Resume" section of [BPS81] one finds: 

... it was found that the inclusion of a dissipative term was much more im- 
portant than the inclusion of the nonlinear term, although the inclusion of 
the nonlinear term was undoubtedly beneficial in describing the observa- 
tions. . . 

Obviously this conclusion is related to dissipation description only and does 
not have the general character. One can find many other evidences in the 
literature which point out the importance of viscous effects. 

Historically, the researchers tried first to include the dissipative effects into 
various long wave models such as Burgers, Korteweg-de Vries and Boussinesq 
equations. There is a vast literature on this subject [Keu48, OS70, KM75, 
Mil76, Mat 76, Kha87, Sug91, Kha97, DD07a, Dut07]. 

In order to include some dissipation into the framework of free-surface po- 
tential flows, we developed the so-called visco-potential formulation [LO04, 
DDZ08, DD07b, Dut07, Dut09]. The kinematic viscosity (or eddy viscosity 
more precisely) appears through local dissipative terms in kinematic and dy- 
namic free-surface boundary conditions. The main peculiarity consists in mod- 
ifying the bottom kinematic condition due to the presence of the boundary 
layer which is assumed to be laminar. Mathematically, this correction proce- 
dure leads to a nonlocal in time term. Note, that from fractional calculus point 
of view this nonlocal term is also a half-order integral. The physical relevance 
of the visco-potential formulation was shown in [LSVO06]. They compared 
model predictions with experiments on the damping and shoaling of solitary 
waves. It was shown that the viscous damping due to the bottom boundary 
layer is well represented by this theory. 

To complete our literature review, recall that there is also an alternative ap- 
proach to potential flows of viscous fluids developed by Daniel Joseph and his 
collaborators [JL94, FJ02, JW04, Jos06]. 

The main goal of the present article is twofold. On one hand, we refine our 
previous dispersion relation analysis. In particular, we do not neglect the evo- 
lution terms it^, where u{t;k) is the wave frequency (5). Consequently, we 
have to deal with complex integro-differential equations. On the other hand, 
we show that bottom boundary layer can induce a disintegrating instability 
of the wave packets. This result is new to author's knowledge. 

The present article is organized as follows. In Section 2 we describe the gov- 
erning equations and perform classical dispersion relation analysis. Then, we 
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derive an equation for the group velocity in Section 2.1. Section 3 contains 
several numerical results and their discussion. Finally, this article is ended by 
outlining main conclusions of the study in Section 4. 



2 Mathematical formulation and dispersion relation analysis 



Consider the 3D fluid domain bounded above by the free-surface z = r]{x,t), 
X = {x, y) and below by the rigid boundary z = —h{x). A Cartesian coordinate 
system with the 2;-axis pointing vertically upwards and the xOy-plane coincid- 
ing with the still-water level. The flow is assumed incompressible and the fluid 
viscous with the kinematic viscosity u. The governing equations of the free- 
surface visco-potential flow have the following form [DD07b, Dut07, Dut09]: 

A(f) = 0, {x,z) efl = R'^ X [-h,ri], (1) 

r]t + Vr]-V(p = (f)^ + 2z/A?7, z = t], (2) 



^ \V(Pf + gV + 2u(l),, = 0, z = r], (3) 



^ / ^-(^^'1^^'"^ dr, z = -h. (4) 
^ 

We can derive corresponding long wave models from the visco-potential for- 
mulation (1) - (4). The derivation can be found in [LO04, DD07b, Dut07]. 

In order to perform the dispersion relation analysis, we have to linearize equa- 
tions (1) - (4) over the flat bottom z = —h. This procedure is classical 
[Sto57, Whi99] and we do not detail it here. Then, we look for the follow- 
ing periodic plane wave solutions: 

(j){x, z, t) = (^(^)e^(^-^-'^(*;^)*) , r]{x, t) = nQe'^k-x~u{tm ^ (5^) 
where k is the wavenumber and u>{t\ k) is the wave frequency. 

Remark 1 It is important to assume from the beginning of the derivation that 
the wave frequency uj explicitly depends on the time t. Consequently, we get 
some additional important terms of the form zt ^^g*'^^ (see equation (6)) which 
were neglected in our previous study [Dut09]. 

We plug the special form of solutions (5) into the linearized version of the 
governing equations (1) - (4). After performing some simple computations 
(details can be found in [Dut07, DutOQ]), we come to the following necessary 
condition of periodic solution (5) existence: 

+ gk tanh(A;/i) - kF{uj, t) {Q'^ tanh(fc/i) + gk) = 0, (6) 
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where we introduced several notations k := \k\, Q := itdtw + iu — 2uk'^ and 
F{uj) is inherited from the nonlocal term: 

^ 

Once the wave frequency 00 is computed from equation (6), the phase speed 
can be immediately deduced by its definition: 

Cp{t] k) := - . 

Remark 2 The classical dispersion relation for water waves u"^ = gk tanh(fc/i) 
can be immediately recovered from (6) if we replace VL by iu and F{uj) by 0. 



2.1 Group velocity 



The group velocity is defined as follows: 

duj(t; k) 



Cg{t;k) :-- 



dk 



There are several physical interpretations of this quantity. One can see it as the 
wave energy propagation speed. Another interpretation consists in viewing it 
as the wavetrain amplitude variation speed. Anyhow, the group speed Cg plays 
an important role in the description of the wavetrains and wavepackets. The 
evolution equation of this quantity can be obtained by simple differentiation 
of (6) with respect to the wavenumber modulus k: 

2nT + gtanh{kh) + gkh{l - tanh^{kh)) - (F + kdkF){n'^ tanh{kh) + gk) 
~ kF ( 2^}T tanh{kh) + {1 - tanh'^ {kh))h + g) =0, (7) 



where fl was defined above and T := itdtCg + icg — Avk. 

Integro-differential equations (6) and (7) are too complex for any mathemat- 
ical analysis. Therefore, we use numerical methods in order to explore some 
properties of the solutions. 



3 Numerical results and discussion 



In order to have an insight into the transient behaviour of the phase and group 
velocities, we have to solve numerically equations (6) and (7) correspondingly. 
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parameter definition value 



V 


eddy viscosity 


10-3 

s 


9 


gravity acceleration 


9.8 ^ 


h 


water depth 


1 m 


At 


time step 


0.05 s 



Table 1 

Values of the parameters used in the phase and group velocities numerical compu- 
tations. 

Hence, in this study we discretize all local terms with a first order implicit 
scheme, while the integral term F{uj,t) is computed in explicit way for the 
sake of computational efficiency. Resulting algebraic equations are solved an- 
alytically. Initial conditions were chosen according to [DDZ08]: 

Fg /I kh \ 

Cr,L r. = \ -r^8Jah(kh) — 2iiyk, cJ. n=( — I r^— ; — r I Re c„|, r, — 4:ii'k. 

The values of all parameters used in numerical computations can be found in 
Table 1 . The time step At is chosen to achieve the convergence up to graphical 
resolution. 



In this section upper and lower images always refer to the real and imaginary 
parts respectively unless special indications are given. Presented here results 
are normalized by ^/gh. Thus, all real parts at the infinitely long wave limit 
kh ^ take the unitary value. 

First of all, on Figures 1 and 2 we show two snapshots of the phase velocity 
at different times. On these Figures we plot also computational results of 
our previous weakly dynamic analysis [Dut09], which consists in neglecting 
differential evolution terms it^ in equation (6). We can constate an excellent 
agreement for the real part of the phase velocity. However, with the simplified 
approach we get only the qualitative behaviour of the imaginary part. This 
discrepancy is located in the region of long waves where the boundary layer 
damping is predominant. The dissipation of short waves is perfectly described 
again. From quantitative point of view, our previous analysis overestimates the 
dissipation rate for long wavelengthes. We refer to [Dut09] for the discussion 
and physical interpretation of the phase velocity behaviour. 



Group velocity snapshots are presented on Figures 3-7. The real part of the 
group velocity Cg is shown only on Figure 7 since its evolution is slow as in 
the phase velocity case. The whole animation of the group velocity can be 
downloaded at [GRP]. 

The imaginary part of the group velocity remains negative until about 1 s (see 
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Figure 1. Real and imaginary parts of the phase velocity at f = 2 s. 
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Figure 2. Real and imaginary parts of the phase velocity at t = 4 s. 

Figure 4). Then, it was very unexpected for the author to observe the forma- 
tion of modes with positive imaginary part (see Figure 5). These destabilizing 
modes may be responsible of wavepackets disintegration. We refer to the next 
section for the discussion of this interesting result. 
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Figure 3. Imaginary part of the group velocity at t = 1 s. 




Figure 5. Imaginary part of the group velocity at t = 2 s. 
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Figure 6. Imaginary part of the group velocity at i = 4 s. 
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Figure 7. Real and imaginary parts of the group velocity at i = 8 s. 
4 Conclusions 



In this work we performed an analysis of linear dispersion relation for visco- 
potential flow formulation. Several numerical snapshots of the phase and group 
velocities were presented. The main particularity in this study is that we 
make no assumption of the weak dependence on time. Therefore, we come 
up with additional terms of the form itu>'{t\ k) which complicate analysis and 
numerical computations. We can conclude that we obtain a good qualitative 
agreement with preceding results [Dut09]. Thus, we validate in some sense 
approximations made anteriorly. However, our computations reveal less oscil- 
latory behaviour in the imaginary part of the phase velocity for intermediate 
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wavelengths. 



Furthemore, we studied the evolution of the group velocity with time for the 
visco-potential formulation. Physically, this quantity represents the energy 
propagation speed. There exists also another interpretation. The group veloc- 
ity Cg{t] k) can be seen as the wavetrain amplitude variation speed. Numerical 
simulations show the appearance of modes with positive imaginary part for 
I < A;/i < I in the group velocity. These modes may destabilize wave pack- 
ets. Physically, this effect comes from the bottom boundary layer. Previous 
numerical study [WLY06] did not reveal this instability since the authors did 
not take into account for the boundary layer effect. Hence, they included only 
a local dissipative term which has a stabilizing effect. 

It is well known that Stokes wavepacket can be disintegrated by the Benjamin- 
Feir instability [Ben67, BF67]. This instability comes from the resonant quar- 
tet wave interaction. The classical result of T. Brooke Benjamin [Ben67], re- 
fined later by R.S. Johnson [Joh77], states that no instability is present if 
kh < ctq and the modes inside a certain band-width grow without bound if 
kh > ctq. The critical value of ctq is estimated about 1.363. There is a certain 
surprising interplay with our results that we would like to point out here. 
Our theory gives less sharper results. That is why we prefer to speak about 
the instability zone defined as U(t) = {kh G M"*^ : lmcg{t; k) > 0} in kh-spa.ce. 
According to our computations, the topological structure of the set U(t) is a 
union of finite number of intervals (see the animation [GRP]). It is curious 
that the critical value ctq hes inside this zone (see Figure 5). 

The wave turbulence can be perfectly described in the framework of potential 
flows. It gives direct cascade of energy from the waves spectral maximum to 
the dissipative region [ZLF92] . However, in the present study we did not inves- 
tigate the effect of viscous interaction with bottom onto the wave turbulence. 
Presumably, it will modify the dissipation region. This question deserves a 
thorough consideration. 

Recently it was shown that the Benjamin- Feir instability can be enhanced by 
dissipation [BD07]. The peculiarity lies in the form of dissipative terms that 
we put into the model. Our computations show the stabilizing effect of local 
dissipative terms [SHC+05, WETGOl, WLY06] and the destabilizing effect of 
our nonlocal dissipation. 

The experiments should confirm our theoretical result, provided that a gen- 
erated wavetrain have sufficiently big wavelength to interact with the bottom 
boundary layer. 
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